clear all;
feature=12533;
classno=11;
data = load('trte11op.txt');
label = data(:,(feature+1):(feature+classno));
data = data(:,1:feature);
[i,j]=find(label==1);
    
        [center,U,obj_fcn] = fcm(data,3);
        %plot(data(:,1), data(:,2),data(:,3),'o');
        %hold on;
        %maxU = max(U);
        % Find the data points with highest grade of membership in cluster 1
        %index1 = find(U(1,:) == maxU);
        % Find the data points with highest grade of membership in cluster 2
        %index2 = find(U(2,:) == maxU);
        %line(data(index1,1),data(index1,2),'marker','*','color','g');
        %line(data(index2,1),data(index2,2),'marker','*','color','r');
        % Plot the cluster centers
        %plot([center([1 2],1)],[center([1 2],2)],'*','color','k')
        %hold off;
